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Abstract. The nontrivial dispersion relation of a periodic medium affects 
both the spectral and the spatial distribution of Cherenkov radiation. We 
present a theory of the spatial distribution of Cherenkov radiation in the far- 
field zone inside arbitrary three- and two-dimensional dielectric media. Simple 
analytical expressions for the far-field are obtained in terms of the Bloch mode 
expansion. Numerical examples of the Cherenkov radiation in a two-dimensional 
photonic crystal is presented. The developed analytical theory demonstrates good 
agreement with numerically rigorous finite-difference time-domain calculations. 
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1. Introduction 

In recent years light-matter interaction in inhomogeneous media, structured on the 
sub-wavelength scale, has attracted a considerable attention. The interplay between 
interference and propagation can result in a nontrivial dispersion relation in such 
a medium. For example, in a periodic dielectric medium [U [2] both dispersion and 
diffraction of electromagnetic waves are substantially modified leading to novel optical 
phenomena, including the emission dynamics modification O [H [5] , ultra-refraction 
[El HIE] an d photon focusing jTUl [Til [12] effects. A nontrivial dispersion relation can 
also substantially modify the Cherenkov radiation |13 [ fT4 l fT5]. For example an electron 
moving in a homogeneous medium with dispersion can emit at any velocity [15] and 
the spatial distribution of the emitted radiation demonstrates intensity oscillations 
behind the Cherenkov cone [TB I fTTl [18] . 

Periodic media modify both spectral [261 [27] and spatial [26] distribution of the 
Cherenkov radiation. At the moment, several studies on the modification of the 
Cherenkov and Smith-Purcell radiation are available. The Smith-Purcell radiation has 
been studied in Refs. |19 | |20 [ l2T j l22 j 123] . The Cherenkov radiation has been used to 
map a band structure of a two-dimentional (2D) photonic crystal in Refs. pUES]. We 
have recently presented a general theory of spectral power modification in both three- 
dimensional (3D) and 2D periodic media [27]. Spatial modifications of the Cherenkov 
radiation have been numerically analyzed in Ref. [26J using the finite-difference time- 
domain (FDTD) method. 

The main purpose of the present work is to develop an analytical theory of the 
Cherenkov radiation in far-field zone in order to provide a simple semi-analytical tool 
to study its spatial distribution pecularities in general 3D and 2D periodic dielectric 
media. To achieve this goal we first derive an analytical expression for the electric 
far-field in terms of the Bloch mode expansion. Secondly we restrict the full k- 
space integration in the Bloch mode expansion of the field to a relatively simple 
surface (contour) integral in the first Brillouin zone of the 3D (2D) periodic medium. 
Thirdly, we restrict this integral further to a simple sum over a small number of Bloch 
eigenmodes. These derived formulas allow us to identify the main contribution to the 
spatial pecularities of the Cherenkov radiation. 

The paper is organized as follows: In Section [2] the general solution of Maxwell's 
equations is summarized. In Section [3] an analytical expression for the Cherenkov 
far-field is derived both for 3D and 2D periodic dielectric media. In Section [4] we 
apply the developed theory to calculate the Cherenkov far-field in the particular case 
of a 2D photonic crystal. Predictions of the analytical theory are substantiated by 
numerically rigorous FDTD calculations. Section [5] concludes the paper. 

2. General solution 

We consider a point (line) charge q uniformly moving with a velocity v in an infinite 
periodic 3D (2D) dielectric medium (Fig. [lj . The medium is described by a periodic 
dielectric function e (r) = e (r + R) with R = J2i ^i a i being a vector of the direct 
Bravais lattice. The U are integers and the at are basis vectors of the lattice. We 
assume also that the medium is linear, nonmagnetic and loss less. The relevant part 
of Maxwell's equations read in SI units: 

V x £(r,i) = - Mo ^K(r,t), (1) 
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r = 0.4a 




Figure 1. A sketch of a periodic medium and charge trajectory. Parameters of 
the infinite 2D photonic crystal used in the numerical example (Sec. HJ are shown. 



V xH(r,t) 



e e{r)^-E{r,t) 



J(r,t), 



(2) 



where E(r, t) is the electric field, H(r, t) is the magnetic field and J(r, t) is the current 
density. In the presence of an arbitrary current density switched on adiabatically 
J(r, t — > — oo) = 0, a general solution of Maxwell's equations JT][2]) for the electric 
field at point r and time t can be given in terms of the Bloch eigenmode expansion 
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d = 2, 3 is the dimensionality of the periodic lattice. The asterisk (*) and ® denote the 
complex conjugate and the outer tensor product, respectively. Ej^J (r) and ^„ (r) 
are generalized transverse and longitudinal Bloch eigenmodes [2, 28J characterized by 



where 



(T) (L) 

the band index n, the wave vector k and the eigenfrequencies u> kn and , 

= 0. The A:-space integration is performed over the first Brillouin zone (BZ) 
of the periodic lattice and the summation is carried out over the different photonic 
bands. Bloch eigenmodes are solutions of the homogeneous wave equation and satisfy 
periodic boundary conditions. These eigenmodes are normalized 

J d d re (r) E%* (r) • E™, (r) = (2vr) d 5 aP 5 m ,5^ (k k') (4) 

and satisfy the completeness relation 



E 

not 



[ d d kV£(r)e(r')E^(r)®E^*(r') 

JBZ 



= (2nf tS id) (r - r') , (5) 

where a, (3 — T or L, and 1 is the unit tensor. 

For a point (line) charge q moving with the velocity v, the current density J{r,t) 
can be expressed as 

J(r,t) =qv8 {d) (r-vt). (6) 
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Using this current density and performing the space integration in the Bloch eigenmode 
expansion |(3|) the following general solution can be obtained 

E(r, t) = - «W £ / d d k f dt> cos Lg? (t - 0" 

(7) 

with ■?) = v/ \v\ being the unit vector in the direction of the charge velocity. Here 
the fact that the longitudinal modes do not contribute to the radiated field [27] has 
been taken into account. In what follows the upper index (T) denoting transverse 
eigenmodes and eigenfrequencies will be dropped. 

3. Cherenkov radiation in the far- field zone 

Expressing the cosine function in Eq. §J§ as a sum of two complex exponential functions 
and using Bloch's theorem Ek n (r) = e kn (r) exp (ife • r), where e kn (v) is a lattice 
periodic function [2], a general solution for the Cherenkov field can be rewritten in 
the form 

v ' n 

where I± is defined as 

I±= f d d k [ dt' e kn (r) (e* kn (vf) ■ v) £ ± (k, t') (9) 

with 

£± (fe, t') = cxp {i [k ■ (r - vt') ± to kn (t - <')]} ■ (10) 

Further taking into account the symmetries of the Bloch eigenmodes, e_fc„ = e kn and 
w_fen = ujkn [2], the following relation for the integrals I± ((9j) holds 

I*_= [ d d k [ dt'e. kn (r)(e*_ kn (vt')-v)£ + (-k,t') 

JBZ J-oo 

= I + . (11) 

Therefore the radiated field |(8]) can be exclusively expressed in terms of the real part 
of the integral /_ 

^r,i) = -^ £>(!_). (12) 

To analyze further the Cherenkov field (fl2|) we limit ourselves to the charge 
trajectories which do not cut dielectric interfaces of the periodic medium. In this 
case "bremsstrahlung" radiation can be neglected and the trajectories themselves 
are necessarily oriented rationally with respect to the periodic lattice. In this case 
e kn i v t') • 'W in J9j) is a one-dimensional (ID) periodic function with a period a defined 
by the orientation of the charge trajectory therefore one can Fourier expand it as 
follows 

e* kn (vt')-v= jr c nm (fe)exp fi^l m A (13) 

in— — oo 
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with Fourier coefficients 



\J d£ (e* kn m ■ v) exp (-i^n*) ■ ( 14 ) 
This allows us to rewrite © in the form 

°o /• rt 

(k)£ nm (k,t r ) (15) 

„JBZ J -00 



m— — cjo 



i ( fc • (r - vt') - uj kn {t - t') + ^^-mt' 



'BZ 

with function £ nm (k,t / ) defined by 

£nm{k, t ) = CXp 

= exp[i/ nm (fc,0]- (16) 

We are interested in the field far away from the trajectory of the charge. In 
the far-field zone the following relation \r — vt'\ A holds for all moments of time 
t' < t. If this condition is fulfilled a small variation of wave vector fc results in 
rapid oscillations of the exponential function £ nm (k, t'). Taking into account that the 
function ekn (r) c nm (fc) is a slow function of the wave vector, the main contribution 
to the integral /_ in the far-field zone comes from the neighborhood of /c-points where 
the variation of the phase / nm (fe, t') is minimal. Such stationary fc-points are defined 
by the relation 

V fe /„ m (fc,t') =0, (17) 
which explicitly reads as 

T — Vt 1 

v g (k) = V fe w„(fc) = — — . (18) 

Where v s (k) is the group velocity of the Bloch eigenmode. Relation lfl8|) can be 
written in the equivalent form 

(t, g (fc)-«).(r-«t) = l^^L | (19) 

with its right hand side being positive for all moments of time t' < t. Taking that into 
account, the integration in 1|15|1 can be restricted to the part of the Brillouin zone, 
BZi, containing all wave vectors whose group velocities fulfill the relation 

(v e (k) -v)-(r- vt) > 0. (20) 

Further, noting that in BZi there does not exist any stationary points for t' > t, 
the ^'-integration in (fl5| can be extended to the whole real axis without severe error. 
With a good accuracy the integral Z_ can be approximated by 

/ d k I dt &kn (^*) Cum 

(k)£ nm (k,t'). (21) 

JBZi J-00 



m=— 00 ' 



Using the integral expression of the Dirac delta function, the ^'-integration can be 
easily performed resulting in 



OO r. 

m=-no J BZi 



(r)c nm (k) exp [i (fc • r - u kn t)} 



2-7T \v\ , 

S uj kn ~k v + L_Lm . (22) 
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0.1 Frequency [a/A] 0.25 

Figure 2. Diagram to illustrate the generalized Cherenkov condition. The 
band structure of the 2D photonic crystal (Fig: [2j for TE (transverse electric) 
polarization as well as the set of planes for m = 0, m = —1 and m = —2 and the 
charge velocity v = 0.15c x are shown. The intersections of the band structure 
with the planes define the integration contour C m (dashed lines). A color coded 
projection of the contour on the first Brillouin zone is also shown. 



Furthermore, using the integral relation 

d d rf(r)S(g(r)) = [ d^V-^, 
Jdv |Vff(r)| 

where dV is a surface denned by the equation g (r) = 0, the d-dimensional integral in 
l(22)l can be reduced to the (d — l)-dimensional integral 



*?1 



exp [i (fe • r - u) kn t)] , (23) 



\v s (k)-v 

which finally gives the following expression for the Cherenkov far-field 



■exp[i(fe-r-o; fen i)])|. (24) 

The integration is performed over the part of the Brillouin zone, BZi, defined by 
(|20p and the integration surface (3D case) or contour (2D case) C m is defined by the 
generalized Cherenkov condition |27) 

ujkn = k-v m \v\ , (25) 

a 

where m is an integer. In what follows we will refer to the integration surface (contour) 
l(25)l as Cherenkov surface (contour). 

Equation ((24)1 is the main result of the present section. We have demonstrated 
that a restriction of the integration range from the whole Brillouin zone to the solutions 
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of the generalized Cherenkov condition ([25| in BZi (|20|l is possible. The generalized 
Cherenkov condition lf25"|) chooses all Bloch modes contributing to the Cherenkov 
radiation. A graphical illustration of the Cherenkov condition (|25)l is presented in 
figure [H The band structure of an infinite square lattice photonic crystal is shown. 
The dielectric constant of the background medium is e — 12 and the radius of the 
air holes is r = 0.4a, where a is the lattice constant (Fig. [J). The manifold of the 
Cherenkov wave vectors (the integration contour C m ) is given by the intersection of 
the band structure, Wfc„, with the set of planes /(fey) = \v\ fen — \v\ ^-m for different 
m. Here fey is the component of the wave vector parallel to the charge velocity. The 
slope of the planes is defined by the charge velocity, being v — 0.15c x in this example. 



3.1. 2D periodic media 

To further simplify the integral (|23|) we can parameterize the contour C m by its arc 
length s. Then the contour integral lf23|) with d — 2 can be transformed into an 
integral over s 

e„ (k(s),r)c nm (k(s)) 



2 '?/ 



ds- 



\v s (k(s))-v\ 

■exp{i[k(a)-r-u n (k(8))t]}. (26) 

The main contribution to the integral comes from the fc-points fe(s") in whose 
neighborhood the phase h(s) — k(s) ■ r — u> n (k(s))t is stationary with respect to 
the variation of s 

d 



k{s) ■ r - u n {k(s))t] 

OS 







F>k u 

(r-vt)-— =0. (27) 



The second equality in lj2"7|) holds as the derivative of the Cherenkov condition ([25)1 
with respect to the arc length s results in 

r)k 

( Vg (k( S ))-v) ■ — = (). (28) 



Combining the stationary phase condition (|27|) with relation l(28|) and taking into 
account the definition of BZi lj20]) . one can see that for the stationary Bloch modes 
the vector v g {k u ) — v must be parallel to the vector r — vt. In other words, only the 
eigenmodes whose group velocities in the coordinate frame moving with the point (line) 
charge pointing towards an observation direction in this coordinate frame, contribute 
to the Cherenkov radiation in the far-field zone. This statement is illustrated in 
figure [U The group velocity contour corresponding to the integration contour C m 
is shown. The color coding is used as in Fig. [21 the same color corresponds to the 
same frequency and consequently to the same wave vector. Main contributions to the 
integral ([26| for two different observation directions are depicted. For the direction Q\ 
only one Bloch mode with group velocity v^ 1 satisfies the stationary phase condition, 
while for the direction # 2 /3 there are two modes with group velocities v^ 2 and v^ 3 
fulfilling the condition. There does not exist Bloch modes satisfying the stationary 
phase condition for forward observation directions characterized by angles smaller than 
01. 
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Figure 3. Group velocity contour corresponding to the integration contour C m . 
The main contribution to the Cherenkov radiation in the direction 6\ ($2/3) comes 
from the Bloch mode(s) with group velocity (velocities) Ug 1 (i)g 2 and Ug 3 ). The 
same color coding as in Fig. [2] is used. 



The integral (|26|) can be further approximated by expanding the phase h (s) near 
the stationary points k(s u ) = fe" in a Taylor series up to quadratic order 



h (s) « k" ■ r - uo n (k v )t + I (r - vt) ■ (s ~ s v f 
and extending the integration range to the whole real axis 

Z_ w 2tt2^ \Zn~TU\ fi exp{i [k v r- uj n {k v )t\\ 



(29) 



\v s (k»)-v\ 



ds exp 



(r — vt) 



d 2 k v 
ds 2 



is- 



(30) 



Where the summation is taken over all stationary points v. The resulting integral can 
be evaluated analytically [29] 



dte ibt = 



which together with the relation 

2uv 



exp 



i^sign(fe) 



(r — vt) 



d 2 k 

ds 2 



sign (r — vt) 



d 2 k v 
ds 2 



\(r-vt)\ 



d 2 k u 



ds 2 



(31) 



(32) 



results in the final expression for the Cherenkov electric field in far-field zone for a 2D 
periodic dielectric medium 
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•expji ^sign ((r - vt) ■ K."} } >. 



(33) 



Here KT 



ds 2 



is the curvature of the contour C m at the stationary point 



k v . Relation (TSl holds because of # _L and (EI 



3.2. 3-D periodic media 

By introducing a 2D coordinate system with unit vectors J^- and tangential to 
the integration surface C m , the surface integral ([23]) with g? = 2 can be expressed as 

„ \^ f f , , dk dk e n (k(s) 7 r)c nm (k(sj) 
J dsi ds 2 \v g (k{s))-v\ 

• exp [i (k(s) ■ r - uj n (k(s))t)} . (34) 

where s G R 2 . Similar to the 2D case, the main contribution to the integral comes from 
the neighborhood of fc-points k(s v ) — k v where the phase h (s) = k(s) ■ r — uj n (k(s))t 
is stationary 



x dk" n 



(35) 



with i € {1,2}. Choosing a local coordinate system at the stationary point k v with 
basis vectors J^- and J^- along the main directions of the surface curvatures the 
following form of the Taylor expansion of the phase h (s) can be used to evaluate the 
integral ([34| 

^ k v ■ r - u n (k v )t 



vt) 



d 2 k v 



(a 



d 2 k v 

~W 2 



(6 - s v 2 y 



(36) 



where and %-h are the main curvatures ¥i\ and K.K of the integration surface 
C m . The integration limit is than extended to the whole real plane. Using twice the 
relation l|3ip the Cherenkov electric field in far-field zone for a 3D periodic dielectric 
medium reads 



E$ D \r,t) 



q\v\ 



2-7r£o |r — vt 

• exp {i [ft" -r-a; n (*:")*]} 

• exp | i 



-sign ( (r - 



I c« (k v , v) c nm (k u ) 

")*]} 



■ exp |i ^sign ((r - vt) ■ } L 



(37) 



with the product of the main curvatures K^KJ^ being the Gaussian curvature of the 
surface C m defined by the generalized Cherenkov condition (|25|) . The vectors ICi are 
defined by & = 
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4. Discussion and numerical examples 

Formulas i|24j) . (j33|) and IpT} constitute the main result of the present work. In the 
far-field zone the electric field generated by a point (line) charge uniformly moving in 
a 3D (2D) periodic medium is dominated by a small number of Bloch eigenmodes of 
the medium. To calculate the far-field, these Bloch modes (their wave vectors) should 
be calculated as a solution of the generalized Cherenkov condition (|25| . In turn, the 
spatial variation of the Cherenkov radiation is dominated (i) by the interference of 
these Bloch modes at the observation point and (ii) by the topology of the dispersion 
relation at the Cherenkov surface (contour) in Eqs. 11241) and l(37|) (Eqs. (|24)) and 

63)). 

In what follows, we apply formulas {24} and (|33| to study the spatial distribution 
of Cherenkov radiation in the 2D photonic crystal depicted in Fig. [TJ The line charge 
oriented perpendicular to the periodicity plane of the crystal moves along the x- 
axis with a velocity v = 0.15c, staying always in the space between air holes. The 
corresponding current density {6} couples only to Bloch eigenmodes with an electric 
field polarized in the periodicity plane (TE polarization). The first TE photonic 
band of the considered crystal is presented in figure [2 The band structure as well as 
group velocities and Bloch eigenmodes were calculated using the plane wave expansion 
method [30] , 

In order to calculate both the electric field lf24|) and its approximation l(33|) the 
set of wave vectors contributing to the far-field should be calculated. This set can be 
found as a numerical solution of the generalized Cherenkov condition (dashed line in 
Fig. [2]). In contrast to the homogeneous medium case, such a solution does exist for 
an arbitrary charge velocity [27] . 

Having the set of wave vectors C m , the Bloch modes efc„ and the Fourier coefficient 
Cnmik) can be calculated. The Fourier coefficients c nm (k) give the coupling strength 
between the current associated with the moving charge and the Bloch eigenmodes. As 
it can be seen from the definition lfl4|) . with increasing index |m| the Fourier coefficients 
become smaller, reducing the contribution of the higher frequencies to the Cherenkov 
radiation in the far-field zone. Finally, the Cherenkov field ([24]) can be calculated by 
direct numerical integration. 

To calculate the stationary phase approximation of the far-field 1)331) , stationary 
wave vectors k v should be calculated for a given observation direction r — vt. This can 
be done by parameterizing the Cherenkov contour by the arc length s and looking for 
all Bloch modes whose group velocities (Fig. [3]) satisfy the stationary phase condition 
(v g (k 1 ') — v) ff (r — vt). Furthermore, the Bloch modes, the Fourier coefficient and 
the curvature of the Cherenkov contour (see Appendix) should be calculated and 
summed for these stationary points only. This reduce the computational demands 
considerably. 

In figure HI a numerical calculations of the stationary phase approximation of the 
far-field is presented for v = 0.15c. The E x (r,t = 0) and E y (r,t = 0) components 
of the electric field are shown in the upper halves of the top and bottom panels, 
respectively. Only contributions from the first three sections of the Cherenkov contour 
corresponding to m — 0, —1, —2 have been analyzed. In the case of a homogeneous 
medium, the standard Cherenkov condition would impose a minimal charge velocity 
above which Cherenkov radiation is possible u m i n > c/n e s ~ 0.457c. Here an 
effective refractive index of the considered periodic medium is asymptotically equal to 
n>efi — y/e e s ~ 2.186. Although the velocity of the line charge is considerably smaller 
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Figure 4. E x (r,t = 0) (top panel) and E y (r,t = 0) (bottom panel) components 
of the Cherenkov electric field for a charge velocity v = 0.15c. The results of 
the stationary phase approximation (1331) and direct numerical integration using 
FDTD methods are shown. 



than Wmin, in the periodic medium the non-evanescent field can be clearly seen far 
apart from the charge trajectory (Fig.|4]). Another characteristic feature of the spatial 
distribution of the Cherenkov radiation for v = 0.15c is a backward-pointing radiation 
cone |26J. The field in the forward direction for observation angles smaller than 9\ is 
zero. 

The zero field within the backward-pointing radiation cone is associated with 
the absence of stationary solutions for observation directions in the cone (Fig. [3]). 
This is a direct consequence of the stationary phase approximation. In figure [5] 
a), the stationary phase approximation of the far-field 10 lattice constants apart 
from the trajectory is compared with the direct numerical integration of the integral 
representation of the field (|24|) . An excellent agreement between these two solutions 
can be seen up to the Cherenkov cone 6*i, where the stationary phase approximation 
breaks down. The cone angle corresponds to the fold in the group velocity contour 
(Fig. [3]). The curvature of the Cherenkov contour is zero at the corresponding wave 
vector and the Taylor expansion (|29| fails to reproduce the contour accurately in the 




Figure 5. Comparison of the integral representation of the far-field Il24|l 
(dashed red line) with the stationary phase approximation l|33|l (black solid line) 
and FDTD calculations (blue solid line). An absolute value of the horizontal 
component of the electric field is shown 10 lattice constants apart from the charge 
trajectory. 



vicinity of this point. Zero curvature leads to the diverging field (Fig. [5]-inset). The 
field calculated with Eq. (|24|) is small but finite in forward directions. 

To substantiate our analytical calculations, a direct numerical integration of 
Maxwell's equations (HE]) using rigorous finite-differences time-domain (FDTD) 
method [31] has been performed. The corresponding FDTD field distributions are 
presented in the lower halves of the top and bottom panels of figure HI 200a x 70a 
lattice of holes with discretization A = a/18 has been used for the FDTD calculation. 
The simulation domain is surrounded by a 3a- wide perfectly matched layer (PML). 
The integration time step is set to 98% of the Courant value. In order to describe a 
continuous movement of the charge, the source ((6]) is modeled by a discrete Gaussian in 
the x-direction with a standard deviation a — A and by an appropriately normalized 
Kronecker delta in the y-direction 

\2" 



j(xiA,yiA) — ^_L fc — exp 



v 



2A 2 

where yu = 35a is the center of the crystal in the vertical direction. In order to 
compare directly FDTD results with the predictions of formulas (|24| and (J33J) , the 
static contribution as well as the higher frequency contributions to the FDTD field have 
been filtered out. An overall good agreement between the results of the analytical and 
direct numerical calculation is obtained (Fig. [4] and E]). The FDTD calculation follows 
nicely the main characteristic of the analytically predicted field. The difference in the 
absolute values between FDTD and analytical calculations (Fig. [5]) can be associated 
with the residual reflections from the perfectly matched layer. 
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5. Conclusion 



In conclusion, we have developed an analytical theory of the Cherenkov radiation in 
the far-field zone. The field far apart from the charge trajectory can be calculated as a 
surface (contour) integral over a small fraction of the first Brillouin Zone in a 3D (2D) 
dielectric medium. We have shown that the main contribution to this integral comes 
from a small and discrete number of fc-points. This opens the possibility to calculate 
the integral approximately, but with high accuracy. We have also shown, that the 
spatial variation of the Cherenkov radiation in the far-field is due to the interference 
of a few Bloch eigenmodes as well as the topological properties of the Cherenkov 
surface (contour). This has been defined as a manifold of all fc-points contributing 
to the Cherenkov radiation for a given charge velocity. Simple formulas have been 
derived for the Cherenkov far-field both in 3D and 2D cases. We have compared the 
developed analytical theory with numerically rigorous FDTD calculations. A good 
agreement between these two methods has been demonstrated. 
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Appendix 



To calculate the curvature using its definition 

d 2 k(s) 



ds 2 



one should numerically evaluate the second derivative of the wave vector k(s) on 
the Cherenkov contour. This is typically associated with a large numerical error. 
In what follows we propose to use an alternative method to calculate the curvature 
which usually results in smaller numerical error and involves calculations of the first 
derivative of the wave vector on the contour and the second derivatives of the dispersion 
relation uj n (k). 

Taking the derivative of (|28| with respect to s we obtain 

d 2 k 



(v s (k(s)) 



dk d 



q s 2 >■ s\ ' ' ds ds 

The second term in IjA.ip can be rewritten as 

dki 
~ds~ 

dk, 



v e (k(s)) = 



(A.l) 



E 



d 9 



E 



u (fc) — 

ds dkjdki ™ ds 
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where Hu, n (k) is a Hessian matrix of the dispersion relation ui n (k). We can now 
rewrite IjA.ip in the form 

d 2 k 
~ds~ 2 



— ■ (v g (k(s)) - v) = M 



or 



K\v g (k(s))-v\ S ign IC-(v g (k( s )) 



= M 



with K. — |y£. Finally the following expression for the curvature is obtained 

K = | J^i r. (A.2) 
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